Exotic Behavior of Heavy-Flavored Meson Matter 
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In this work, we study the thermodynamic behavior of heavy-flavored meson matter in the frame- 
work of (a, w)-meson-exchange model in relativistic mean field theory. We find a decreasing of the 
effective masses of D and B mesons as the temperature increases. By using the effective mass and 
maximum value of dissociation temperatures available from lattice QCD, the masses of the bound 
states DD and BB are estimated in 2 MeV for both molecules. For the B-meson matter, the pres- 
sure presents an exotic behavior, being negative for temperatures above 6.6 times the deconfinement 
transition temperature T c . In addition, the ratio of pressure to energy density is similar to the value 
predicted for systems that behave as dark energy matter. 



I. INTRODUCTION 
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In 2003 the Belle Collaboration [l| reported the first experimental evidence for non-qqmesons. Thence- 
forth several new states (called X, Y and Z) have been discovered in B-meson decays with masses 
between 3.9 GeV and 4.7 GeV for charmonium sector and Zf,(10610) and -Z&(10650) for bottomonium 
sector. These states, however, do not fit as excited states of charmonium described in the framework of 
constituent quark model [6]. Furthermore, they have an exotic decay modes, in which charmed mesons 
are suppressed. 

Thus, a large amount of effort has been directed to understand the structure of these exotic mesons. 
In particular, one new meson that has strong evidence of its existence is the V(3872), which has recently 
been measured in LHCb 5]. The quantum numbers that have been sugg ested by PDG Q are isospin 
1 = and J PC = 1 ++ or J = 2 h . From a theoretical point of view, several proposals have been 
■ attempted in order to provide some insight in this issue, as models based on glueball [8, 9], mixed state 

(cc — DD*) [l0|, tetraquark state [TO] and molecular state (D°D*° + c.c.) [13]. 

Notice that the notion of hadronic molecular states is not new, being proposed several decades ago 
fl3l [HI- Nevertheless, it has been largely employed to get some understanding about the properties 
of exotic states, mainly when the observed meson masses are rather close to a respective two-meson 
threshold. In the example discussed above, the mass of X(3872) is very close to DZ)*-state threshold. 
The same idea is used to other exotic states, with possible molecules like D*D*, BB* , B*B* and so on. 

Pursuing in the molecule approach, another interesting possible combination of heavy-flavor mesons is 
a bound state between D and D, called X(3700). This state has spin-0 and isospin-0 and was predicted 
by Gamermann et al. [HJ], but still here there is no experimental evidence for this state. However, BES 
Collaboration [T(| [TtJ observed in the process e + e~ — > DD two states: the "expected" ^(3770) and 
a probable new state in the energy region between 3.700 and 3.872 GeV. In addition the results from 
branching fraction for 0(3770) — > DD cannot be explained if "0(3770) is a pure cc state, but if it contains 
four-quark admixture [HI, ■ 

Many works have been devoted to study of bound state of D and D that support the existence of 
X(3700) [ll[20|-[23L In these analyses different frameworks were used, as the couple channels [lH, Bethe 
Salpeter equation |20j| , meson exchange model [2l[ , chiral quark model [22[ , QCD sum rules [23| and mean 
field approach [24j |. However, it is worth remarking that there is in literature divergent understanding of 
X (3700) as a molecular state [25|. Also, it deserves mention the meson exchange approach discussed in 
Ref. [2lJ, in wh ich is argued that the vector meson exchange plays an important role to create S wave 
bound states. On the other hand, in canonical quantization at mean field approximation [2~ij the dynamic 
solution for the vector field is zero; the effective mass of D-mesons is constant for low temperatures and 
almost zero only at very high temperature of T = 1.2 GeV, in which the D mesons should be dissociated. 
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The usual deconfinement transition temperature for two flavors dynamical quarks is T c — 172 MeV 
[26j , on the other hand lattice QCD in the quenched app roximation is estimated the value of the decon- 
fmement transition temperature which is 262.5 MeV [23]. In terms of the usual value of T c the suggested 
dissociation temperature in Ref. [HI is then T — 1.2 GeV s» 6.7 T c . Calculations in lattice QCD, however, 
suggest that temperature of charmonium dissociation is a value between 1.6 T c and 2.35 T c [28|, [29[. In 
addition, a recent QCD sum rule result [30j predicts at temperature between 1.52 T c and 1.68 T c the 
charmonium should be dissociate. 

The effective mass of the hadron is the mass in which it seems to carry in the medium. If the effective 
mass is lower than rest mass, the hadron is found as bound state among it and the medium. The 
interaction of nucleons in nuclear matter have been studied using the Walecka model [31(. At finite 
temperature the interaction of many body nucleon-anti-nucleon system presents the effective mass of the 
nucleon is almost zero at 1.74 T c and a sea ofpairs of N-iV are produced in the nuclear medium and this 
system behaves as massless bosons system [32|. Another many-body cases are studied in the interaction 
of mesons with nucleus, such as bound state among meson and nucleus are called mesic nuclei has been 
explored in different ways: D-nucleus [33-35] J/^-nucleus [IB] p-nucleus [13], K- nucleus [Hj]. In these 
situations, the effective mass also appears as a relevant parameter to study the bound states, and has 
been used in Refs. [H, |35T - |37| , where if effective meson mass decreases in the nucleus so mesic nuclei 
exist. 

Hence, taking as motivation the analysis of exotic states within hadronic molecular approach and meson 
exchange models discussed above, in this work we apply the relativistic mean field theory developed in 
Refs. [24], [3l|, [HJ to study the thermodynamic behavior of a system constituted of heavy- flavored mesons 
in the framework of path-integral formalism. We use the following idea: whether the interaction among 
mesons builds a simple structure likes a di-hadron molecule, as suggested by Refs. [TH [20l - [23l ] . so a many- 
body system of hadrons at high temperatures can form more complicate structures like a multi-meson 
system. In special, we analyze the effective mass of the D meson, m e ff, where if m e // decreases we have 
a bound of D in D mesons matter or should have some bound states of D — D or X(3700) in meson 
medium. In addition, if ra e // goes to zero, a sea of pairs of DD or a multi-meson strongly interacting 
system is formed. 

The organization of this paper is as follow. In Sec. II, we introduce the Lagrangian density that describes 
a system of D-mesons and calculate the relevant quantities to study its thermodynamic behavior. The 
analysis of thermodynamics of charmed meson matter is discussed in Sec. Ill, while in Sec. IV this 
approach is extended to bottomed meson matter. Finally, Section V presents some concluding remarks. 



II. THE FORMALISM 



We start by introducing the effective Lagrangian density that describes a system of D-mesons inter- 
acting with scalar (a) and vector (u>) mesons [39(, 

C = (d^D)(d"D^)-m 2 D DD^ - -F, W F^ + ^m^w" 

+\{d ll &){d»a) - l -mla 2 +g DDa DD^ + ig DD ^[Dd^ - (d^D)D% (1) 

where F^ v — d^u) v — d^ut^, D = (D°, D + ) is the doublet of the D meson field, nip, m w and m a are the 
respective masses of D, a and u> mesons, and g^Do an d 9ddu are t ne coupling constants for the DDa 
and DDu interactions. 

In Lorentz gauge, i.e. O^uj^ = 0, the equations of motion obtained from Eq. ([IJ are 

d^a + ml* = g DDa DD\ (2) 
d^^ + mlu" = -ig DBuJ [Dd v rf -(d"D)D% (3) 
d'^D + (m 2 D -g DDa a + g^Lo^^D = 0, (4) 

where 8^ = 3^ + ig D r>^^ 

For our purposes, it is convenient to work in the context of the mean- field approximation [24j, |3l|, 
replacing a and ui meson fields by classical fields, 



a = {a 



<-°) , (5) 
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with w M = for p, ^ 0. This approximation allows us to rewrite the Lagrangian density in Eq. ((T|) as 
C = D [-d^ - {ml - g DDa a)} £>t + ^m 2 w 2 - l m 2 a a 2 

-9ddu & 3°, (6) 

where j° = i[Dd°D' i — (d°D)D>], Besides, it can be seen from Eqs. ©-® in mean-field approximation 
that a oc p s and Q oc p, with p s and p being the scalar and number densities, respectively. 

To study the thermodynamic behavior of the system introduced above, we assume that it is in ther- 
modynamic equilibrium and at finite temperature T and density p, of Z?-mesons. In this sense, we use 
imaginary time (Matsubara) formalism, and the grand partition function takes the form [40j . 



Z oc 



J VD VD 1 cxp j- jf 9 dT J d 3 x [C 



E + Pjo] 



(7) 



where j3 — 1/T and p are the inverse of temperature and chemical potential, respectively, and Ce is the 
Lagrangian density given by Eq. ^ in Euclidean space [50j . Then, after the integration over the fields 
D and D' the obtained thermodynamic potential becomes, 



fi(T,/i) = V 



1 2-2 1 2-2 

—m„a — —to, ,U) 



V 00 

-- y 



d 3 p 
(2^)3 



In 



27m 



iPeff 



+ p 2 +m 2 eff \, 



(8) 



where the effective mass and effective potential of the D mesons in hadronic medium read, respectively, 

(9) 



m lff = m D~9DDaV, 



Veff 



V is the volume. 

In order to obtain the thermodynamic potential in a more tractable form, we use zeta-function regu- 
larization techniques [40M43| . In this scenario, Eq. (JS| can be rewritten as 



orf./O V ( \m 2 a o 2 - X -mlG? 



where the zeta function £(s) is given by 



(10) 



COO 



E 



r s- 



d 3 p 
(2Vp 



/ 2wn 



Weff +P + m eff 



(47r) 3 / 2 r (s) 



E 



n— — oo 



27m 



iPeff 



l eff 



Notice that £ '(s) denotes the derivative of £ with respect to the argument s. 
Then, by using the identity 



+oo _ I — 



a T(r/) 
4tt 1 ' 



/ar(7j) 



cos (2irnc) 



'a q 



K„ 



{ 2imq\ 



and after some manipulations the thermodynamic potential in Eq. (|10[) becomes 



n(T, p) = V I im CT a 2 - 1 



— > cosh(/3n^ 

71- Z { J 



eff) 



n=l 



™eff 

n/3 



K 2 {n(3m e ff), 



(11) 



(12) 



(13) 
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where K v is the modified Bessel function of second kind, and U vac is the T, //-independent vacuum 
contribution which comes from the pole term in Eq. (1121) . 



1 



64tt 



{2r 



l eff 



A 



A 4 + 2mi ff In (m^) + 2 (A 4 - m^) In (m^ + A 2 ) } 



(14) 



with A being a cutoff. We will omit this term henceforth. 

Another requirement necessary to the thermodynamics of this model is the analysis of gap equations, 



= 

da 

dn - 
dw 



(15) 
(16) 



These expressions give the values for a and u> that extremize the thermodynamic potential 51. So, the 
use of Eq. (JT3J) in gap equations (fT5l) and (p~6|) yields 



9dd<j 

~P S 

9dDuj 



hi 



where the scalar and number densities are written as 



' 71=1 

m 2 00 



^ cosh {Pn^ef^-Kx (n/3m e ff ) , 



P 



Mr X smh(f3nfi eff )^K 2 (nl3m eff ). 



(17) 
(18) 

(19) 
(20) 



As remarked before, we have omitted the contribution coming from vacuum term. 

It is worthy mentioning that if the numbers of D and D mesons are equal in hadronic matter, which 
means a vanishing effective chemical potential "e// = 0, from Eqs. (fT5|) and (|19p we see that p = and 
therefore Q = 0, just as the situation reported in Ref. [241 ] . 

Now we derive other relevant thermodynamic quantities from the thermodynamic potential given by 
Eq. (JT3J). For example, the pressure reads 



an 

dV 



1 -2 1 -2 

= ^rn^u) — —m„a — U v 



^ oo 

— ^cosh (J3np eff ) 



n/3 



K 2 (nf3m ef f), 



(21) 



Therefore, the values of a and u that extremize the pressure must also satisfy Eqs. (IT71) and {18} . 
In addition, the entropy and energy densities at chemical equilibrium are given by, 



s(T) 



dp 
df 



^2 X] 



3 



"eff 



K 3 (n/3m e ff), 



(22) 



and 



e(T) = (- P + Ts)^ ff=0 

1 _o 



^ OO 



n/3 



K 2 {n(3m eff ) + — ^ifi(n/3m e //) 



(23) 



respectively. 
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III. EFFECTIVE MASS OF D-MESONS 



We now analyze the thermodynamic behavior of the model introduced in previous section under the 
change of values of the relevant parameters. We consider the situation fi e ff — 0, where the number of D 
and D mesons are equal in mesonic matter. We use the following parameters: = 1.87 GeV, m a = 0.5 
GeV 7], g DDa = 2.85 GeV Q and deconfinement temperature T c = 172 MeV [26lj. 

In Fig. Q] is plotted the sigma field in Eq. (fl7|) as a function of temperature. We note that a does not 
exist for temperatures above of T = 1.34 GeV. 

Besides, a has a maximum value at T = To ~ 1.15 GeV or To rj 6.68T c , which is the critical value in 
which we have an almost vanishing effective mass. This critical temperature can be obtained analytically 
by taking m e ff ~ in Eq. {jTTJ) and using the asymptotic formula of Bessel functions, that is 



Tr 



m D m a 
9dd<j 



(24) 



Replacing the values of the parameters above mentioned, we get Tjj 
with the Fig. [Q 



1.137 GeV or 6.6 T c , in agreement 



D 

O 

t> 0.6 




0.0 0. 



FIG. 1: Plot of a in Eq. (|17[l . as a function of temperature, at chemical equilibrium. 



Let us now turn our attention to the effective mass of the D mesons in hadronic matter. We plot in 
Fig. [2]_the values for m e ff as a function of temperature which we solve the gap equation (fl7|). As in 
Ref. [24]], we can see that the effective mass reduces as the temperature increases, and a transition from 
interacting gas to very strongly interacting matter appears at Tjj. Moreover, this temperature is too high 
compared with the temperature of charmonium dissociation, where its temperature is a value between 
1.6 T c and 2.35 T c [28|, So, in this scenairo the £>-mesons should be dissociated. 

We can estimate binding energy of a meson molecule composed of two mesons D and D, defined by 

e(T) = m DD -2m eff (T), (25) 

where m DD = 2 mo and 2m e //(T) is the mass of £>Z)-bound state. Then, considering the maximum 
temperature of 2.09 T c in Fig. [2j the effective mass has a value 1.869 GeV. Therefore, a DD-meson 
molecule has a bind energy of 2 MeV. For the sake of comparison, notice that in Ref. [2l| is reported the 
binding energy of 1.4 MeV for the S wave DD bound state at T = (with an cutoff A= 1.5 GeV and 
m a = 400 MeV). Thus, this result is close with our findings, but in our formalism the bound state only 
could be appear at T > 0. Note also that if we consider m a = 400 MeV at temperature 1.92 T c we get 
the same effective mass of 1.869 GeV. 

In Fig. ([2]) is plotted the effective mass as function of temperature by taking bigger values of the 
coupling constant g^Qa with respect to the vacuum value of g^Da- ft can be remarked that the critical 
temperature decreases as goQa increases, as explicited in Eq. (f2"4"]l . For instance, with g DDr7 = 9 GeV 
in Eq. {22) we get T D 0.36 GeV or 2.09 T c , in agreement with Fig.©. This fact explicits the role of 



G 



the sigma field in binding the D mesons, and the increasing of the magnitude of interaction among the 
sigma and D-meson fields forces the system to experiment a phase transition at smaller temperatures. 
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FIG. 2: Plot of effective mass of D mesons in Eq. (|17(l as a function of temperature, at chemical equilibrium, for 
different values of coupling constant. 

For a better understanding of the behavior of this system, in Figs. [3] and [4] are plotted the thermo- 
dynamic potential written in Eq. (fT5|) as function of effective mass for g D jj a — 2.85 GeV and 9 GeV, 
respectively. As suggested by the previous figures, we observe a first-order phase transition occurring at 
T w 1.31 GeV for g DDa = 2.85 GeV, and T « 0.582 GeV for g DDa = 9 GeV. 




-0.002 - 



0.0 0.5 1.0 1.5 

m e ff (GeV) 

FIG. 3: Plot of thermodynamic potential of D mesons in Eq. (|13|) as a function of temperature, at chemical 
equilibrium and for g D ua = 2.85 GeV. Dashed, solid and dotted lines represent the cases for T = 1.300, 1.313 
and 1.320 GeV, respectively. 

In Fig.[5]is plotted the pressure in Eq. (IT71) as a function of temperature for g D Q a — 9 GeV. We note an 
exotic behavior for values T > 2.09 T c , in which the pressure will assume negative values. On the other 
hand, as remarked in Ref. [24]], for g^DG = 2.85 GeV the pressure has only positive values. Thus, this 
result suggests that a greater interaction between the D-mesons induces a phase transition from gas into 
a highly interacting matter with negative pressure. Another interesting fact is that the ratio of pressure 
to energy density given in Fig. [5] is in a range of values typical to that predicted for matter that has a 
behavior similar to dark energy, i.e. —1 < 2 < —1/3 [44 145| . 
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m efs (GeV) 

FIG. 4: Plot of thermodynamic potential of D mesons in Eq. (| 13f) as a function of temperature, at chemical 
equilibrium and for g D £,cr — 9 GeV. Dashed, solid and dotted lines represent the cases for T = 0.575, 0.582 and 
0.588 GeV, respectively. 



0,010 



-3 



-0,005 




FIG. 5: Pressure of D-meson matter as a function of temperature at gup^ = 9 GeV. 



In addition, another exotic feature of the system with g^Da = 9 GeV is that the energy per pair of 
mesons at temperature To has a phase with an energy of 3 GeV less than the gas phase, as can be seen 
from Fig. [7] 



0,1 0,2 0,3 0,4 0,5 0,6 

T(GeV) 



FIG. 6: Ratio of pressure/energy density of D-meson matter as a function of temperature at g D Q a = 9 GeV. 
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FIG. 7: Energy per pair of D-meson matter as a function of temperature at g D jjcr = 9 GeV. 



IV. EFFECTIVE MASS OF J3-MESONS 

Taking as a guiding principle the heavy quark symmetry, a natural question that appears is about the 
extension of the approach discussed in previous Sections to systems constituted of other heavy-flavored 
mesons. Then, in the following we study the bottomonium analogously to the charmonium situation 
analyzed above. To do so, we just replace the D-mesons with the corresponding B-mesons in the model 
displayed in Eq.([T|): the doublet of J5-mesons is B = (B~ B°). In the present case we use the following 
parameters: m B = 5.279 GeV Q and g BS(T = 8.04 GeV [39|. 

Thus, using Eq. (|24| adapted to i?-meson system, we get Tb ~ 1.137 GeV or 6.6 T c , i.e. equal to the 
temperature of the D-mesons To- It deserves mention that this is so because the coupling constants 
9dd<t an d 9bb<j we have used are proportional to the mass of the respective heavy meson, as shown in 
Ref . [39| , which leave the temperatures Trj and Tb independent of the mass of heavy- flavored meson. 

Furthermore, we note that the temperature Tb is much higher than the values suggested in literature 
for the dissociation of bottomonium, which are: 2.06 T c [46[ and 4.18 T c |47H48j. 

In Fig. [5] is plotted the effective mass as a function of temperature. We observe the effective mass of 
B meson goes to zero at temperatures close to the value Tb- If we consider the maximum temperature 



9 



given by 3.953 T c in Fig. [5]thc effective mass has a vaiue of 5.278 GeV. Here, a meson molecule composed 
of two mesons B and B has a binding energy given by 2(ms — m e ff) or 2 MeV, which is equal to that 
obtained for the D£)-molecule. On the other hand, it does not agree to the results of Ref. [2l| which 
provide a value of 60.7 MeV for the bound state BB. In our case to obtain a value of 60 MeV for this 
bound state, we must consider that the molecule survives at temperature 6.16 T c , which is slightly lower 
than the critical temperature Tb- 




FIG. 8: Plot of effective mass of B-meson as a function of temperature, at chemical equilibrium, for different 
values of coupling constant. 

In Fig. [9] is plotted the pressure as a function of temperature. We note that the pressure has the 
negative values for T > Tb, even working with the coupling obtained in vacuum. This result suggests 
that the interaction among the B-mesons at temperature T > Tb is stronger than the situation at 
zero-temperature. 




0,8 1,0 1,2 
T(GcV) 



FIG. 9: Pressure of the system of B-mesons as a function of temperature at g B g a = 8.04 GeV. 



The ratio of pressure to energy density is shown in Fig. [TO] Differently to the case of D-mesons, with 
the vacuum coupling constant (gsBa = 8.04 GeV) this ratio present values in the range predicted for 
matter that behaves as dark energy one. 

Finally, another exotic feature can be pointed from Fig. II H in which is plotted the energy per pair of 
mesons versus temperature at vacuum coupling constant. We see that the energy per pair of mesons at 
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FIG. 10: Ratio of pressure/energy density of B-mesons as a function of temperature at g B g a = 8.04 GeV. 



temperature Tg has a phase with an energy of 8 GeV less than the gas phase. 




FIG. 11: Energy per pair of B-mesons as a function of temperature at gBBa = 8.04 GeV 
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V. CONCLUDING REMARKS 

In this work the thermodynamic behavior of a system constituted of heavy-flavored mesons has been 
evaluated, with the DD or BB interaction being studied within a a, w-meson exchange model. At mean- 
field approximation, we have shown that the omega meson does not contribute when the system has 
the same number of heavy mesons and anti-mesons. On the other hand, the sigma field has nontrivial 
solution coming from gap equation, and grows as the temperature increases. 

Using the value of coupling constant typical of the interaction between D and D in vacuum, g^Bo- = 2-85 
GeV, our results suggest that the effective mass decreases as temperature is raised. Also, a bound state 
with energy of 2 MeV can be formed at temperature of 2.09 times the deconfinement temperature T c , 
which is slightly above the values predicted in literature for the charmonium dissociation. 

In addition, an interesting result has been obtained for bigger values of vacuum coupling constant, 
where above the temperature of 2.09 T c , we have obtained a strongly interacting phase having a negative 
pressure, and a ratio of pressure to energy density in the range of values typical of systems that present 
behaviors like dark energy matter. Notice, however, that since the coupling constant g£,Q a is not free, 
we believe that the charmed mesonic matter should not show negative pressure. Nevertheless, this result 
might be useful in models of dark energy in which interacting fields could be treated in usual way, as an 
alternative for instance to the noncommutative field theory [49j. 

We have also analyzed bottomed mesonic matter analogously to the charmed system. With the value of 
coupling constant gBBa m vacuum, at temperature of 3.953 T c the bound state BB has a binding energy 
of 2 MeV. This result is lower than the result of Ref. [2l|, which provide a value of 60.7 MeV. For a value 
of 60 MeV for the bound state in our case, we must consider a molecule survives at temperature 6.16 T c , 
which is above the temperature of bottomonium dissociation available in literature: 4.18 T c (47l. [iij. 

For the case of matter of B-mesons at temperature up to 6.6 T c , we again obtain a phase with negative 
pressure and ratio of pressure to energy density within the range of values expected for systems with 
behavior like dark energy matter. It is relevant to notice, though, that the negative pressure can be 
discarded if the bottomonium dissolve at lower temperatures than 6.6 T c , as have been predicted by 
Refs. dfjHU. 
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